Understanding scale invariance in a minimal model of 
complex relaxation phenomena 

P.I. Hurtado^, J. Marro and P.L. Garrido 
Institute Carlos I for Theoretical and Computational Physics, and 
Departamento de Electromagnetismo y Fisica de la Materia, 
Universidad de Granada, E-18071-Granada, Espana 

February 2, 2008 



Abstract 

We report on the computer study of a lattice system 
that relaxes from a metastable state. Under appropri- 
ate nonequilibrium randomness, relaxation occurs by 
avalanches, i.e., the model evolution is discontinuous 
and displays many scales in a way that closely resem- 
bles the relaxation in a large number of complex sys- 
tems in nature. Such apparent scale invariance simply 
results in the model from summing over many expo- 
nential relaxations, each with a scale which is deter- 
mined by the curvature of the domain wall at which 
the avalanche originates. The claim that scale invari- 
ance in a nonequilibrium setting is to be associated with 
criticality is therefore not supported. Some hints that 
may help in checking this experimentally are discussed. 

Key words: nonequilibrium, relaxation, metastable 
states, avalanches, scale invariance, 1// noise. 

1 Introduction 

Scaling behavior described by a simple power law 
is ubiquitous in nature. This is endowed a great 
theoretical interest on the assumption that it re- 
veals prevalence of some underlying feature 
Power-law distributed events have no typical size, 
known as "scale invariance" . The spectral density 
also reveals multiplicity of scales in many time se- 
ries, i.e., the mean square fluctuation goes inversely 
with frequency, known as "1// noise". 

The concept of scale invariance originated in 
equilibrium statistical mechanics. This predicts 
dramatic extension of correlations near critical 



points; any microscopic spontaneous fluctuation 
then triggers events of any possible size with the 
same cause. Thermodynamic equilibrium is a spe- 
cial, pathological case, however, and deep under- 
standing of why complex — out of equilibrium — 
systems are capable of events of any size, e.g., 
crashes in the stock market or disastrous earth- 
quakes, is lacking. 

Most of the cases already studied in detail, 
mainly those in the realm of physics, has now 
a model explanation. For example, the struc- 
tured noise perceived in amplifiers besides a ran- 
dom spectrum jHj is explained based on the physics 
of the electron transport in a vacuum tube [Hj; 
and the Barkhausen noise ^Uj is associated with 
the discontinuous motion of domain walls between 
pinned configurations in a disordered medium 
Specific models do not explain, however, 
the observed ubiquity and universality of the phe- 
nomenon. Among the interesting approaches that 
investigate this aspect, we mention SOC or self- 
adaptation into a critical condition ; the possibil- 
ity that a system naturally lies on the edge between 
order and chaos |12j ; proximity to an standard crit- 
ical point 1 181 1141 ITS] ; and the hypothesis that nat- 
ural selection induces evolution towards a "highly 
structured state" JU- 1 There is no full agreement 
yet on a globally coherent explanation, however. 2 

We present here a new effort towards better un- 
derstanding this problem. We analyzed in detail 
a minimal model of relaxation phenomena which 



iMore specific mechanisms have also been proposed, in- 
cluding those in Refs. lT7) -|25 | . 

2 See the criticism in, e.g., Ref s. ISI 131 IST1 1351 1251 IS7I 1251 . 
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shows not only many scales but also some of the 
basic processes that typically characterize the nat- 
ural phenomena of interest. In particular, under 
appropriate conditions, the model evolution pro- 
ceeds via avalanches, each following a time plateau, 
and this occurs in a way that closely resembles the 
discontinuous variations with time of certain sig- 
nals. For example, a similar relaxation process has 
been described concerning the current through re- 
sistors |21) . the magnetization while a varying field 
induces domain rearrangements (23 HEBE], the en- 
ergy released in earthquakes PIE], and the erosion 
of rocky coasts ■ 

The versatility of the model allowed us to deter- 
mine the conditions in which avalanches are power- 
law distributed. Analysis of these cases reveals that 
the decay (from an initial metastable state) rather 
consists of successive exponential relaxations, each 
with a characteristic well-defined size. We demon- 
strate that a broad distribution of sizes occurs, but 
that this is not to be associated to long-ranged cor- 
relations but to certain randomness. This provides 
some microscopic support to the old suspicion that 
observed electronic and magnetic 1// noises simply 
consist of a superposition of many different typical 
scales, each with a different origin. This was argued 

in Refs.QTIdimEn]; S6e als ° ^ &nd references 
therein. 

The fact is that our model describes a random 
combination of event sizes which produces an ef- 
fective situation which is reminiscent of (equilib- 
rium) critical behavior. Our study does not sup- 
port, however, a critical condition — neither chaotic 
behavior — as a source for the many scales in the 
model. This may have important practical impli- 
cations. 

In spite of its mathematical simplicity, our ex- 
planation for the observed many scales is not phys- 
ically trivial. That is, we describe a time plateau 
before each avalanche takes place which is to be as- 
sociated to entropic metastability. On the other 
hand, we conclude on the existence of dynamic 
(non-critical) correlations. Furthermore, it ensues 
that the many scales definitely require randomness 
in a nonequilibrium setting. 

We claim that our findings are consistent with 
some recent observations, and that it should be 
possible to test them in purposely designed exper- 
iments. 



2 Model and its motivation 

Let the square two-dimensional lattice with binary 
spin variables at the N sites inside a circle of radius 
R. Interactions are according to the Ising ferromag- 
netic energy function, H = - J2{ij) s i s j ~ ^J2i s i> 
where the first sum is over nearest-neighbor pairs 
of spins. Assume open circular boundary condi- 
tions, i.e., bonds leaving the circle are broken. Any 
configuration s = {si — ±1} evolves stochastically 
with time by spin flips with rate: 

g-Affi/T 

w (si -> -Si) = p + (1 - p) — e _ AHi/T , (1) 

where AHi stands for the flip energy cost (we set 
Boltzmann constant to unity). This amounts to 
perturb at random, with probability p, a canoni- 
cal tendency to the equilibrium state correspond- 
ing to temperature T and energy H. This is an 
efficient, non-trivial way of implementing a com- 
plex nonequilibrium situation. That is, dynamics 
involves a conflict — between finite and "infinite" 
temperature — which asymptotically drives the sys- 
tem for any p > towards a nonequilibrium steady 
state. This essentially differs from the equilibrium 
situation for p = |33j . Consequently, similar mod- 
els with competing temperatures have been studied 
during the last decade as a paradigm of systems far 
from equilibrium |34| . 

The singular behavior that ensues for p =/= (in 
the presence of open boundaries) is, in fact, our 
main motivation for studying this system. That is, 
setting p > 0, even as small as p ~ 10~ 6 , induces 
a series of successive time plateaus or short-lived 
"halt states" during the evolution that are not ob- 
served for p = 0. Therefore, £Q) in practice provides 
one of the simplest scenarios one may think of for 
analyzing in detail a whole class of relaxation phe- 
nomena. That is, a principal feature of both this 
simple model and the cases mentioned in section^ 
is that relaxation proceeds by nonequilibrium vari- 
ations that are realized as jumps between locally- 
stable states. 

The present system may be interpreted — though 
this is not essential to the conclusions below — as 
an oversimplified model of a small ferromagnetic 
particle. This, which is relevant to the technol- 
ogy of dense magnetic media [33] 0H1 EH > requires 
one to deal with a large surface/ volume ratio and, 
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consequently, with impurities. In actual speci- 
mens, these typically cause perturbations, e.g., dif- 
fusion of defects will dynamically disturb the lo- 
cal fields, which one may ideally represent by the 
random term in J}. In fact, a similar ansatz 
has already help the understanding of ionic dif- 
fusion during magnetic ordering and it has 
been useful to model microscopic quantum tunnel- 
ing |39| and non-localized perturbing interactions 
and fields 

Whichever the specific interpretation is, the 
model contains a microscopic random perturbation 
which drives it out of equilibrium. This nonequi- 
librium randomness, which is likely to characterize 
also the phenomena of interest (fault slips, elec- 
tron transport, magnetic domain arrangements,...) 
happens to induce the most interesting behavior 
during relaxation. In order to show this, we per- 
formed a series of computer simulations. They 
typically begin with all spins up, Sj = +1 Vi at 
t = 0. For a negative value of the field h, this is a 
metastable state. In fact, for a low enough value of 
the temperature T, the stable state corresponds to 
to = N^ 1 J2i s i — ~ 1- That is, most of the spins 
need to flip to point down along the field direction 
in any stable configuration. For appropriate values 
of the parameter set (T, h, p) - see below, such a 
decay consists of a sequence of transitions through 
short-lived states. This is illustrated in figure 1 for 
a particle of approximately 10 3 spins at T = 0.25 
(which is about 1/10 the critical temperature of the 
corresponding infinite system), h = —0.1 (1/10 the 
exchange energy), and a very small value for the 
perturbation, p = 10 -6 . 

The same behavior ensues for a broad range of 
values for (T, h, p). To be more precise, one needs 
that both domain walls and clusters are well de- 
fined. Otherwise, the jumps are difficult to be ob- 
served and/or one obtains poor statistics. In order 
to assure compact configurations and clusters, it 
turns out necessary a sufficiently low choice for T. 
On the other hand, the parameter p can take a con- 
siderable range of values, provided that its effects 
are comparable to the ones from other stochastic 
sources. This was observed to occur already for 
p = 10~ 6 , which is the value we used in most of 
our simulation runs (see for other choices). The 
other model parameter, h, just aims at producing 
metastability, so that only its sign is really relevant. 
Summing up, the behavior described here is robust 




Figure 1: Typical evolution in which the magneti- 
zation is observed to decay by jumps to the final 
stable state. This is for a single particle of radius 
R = 30 (~ 10 3 spins) at low temperature, and for 
small values of h and p (see the main text). The 
time axis shows t — To in MCS (Monte Carlo steps 
per site) with To = 10 30 MCS; this is of order of the 
duration of the initial metastable state. The inset 
shows a significant detail of the relaxation. 



within a broad region of parameter space, so that 
no fine-tuning of parameters is needed. 



3 Some details of relaxation 

Taking the MCS as the relevant, macroscopic time 
scale, one identifies (e.g., in fig.l) strictly mono- 
tonic changes of m (t) that we call "avalanches" . 
Let us define the avalanche duration A t = \t a — tb\ 
and size A m = \m(t a ) — m(tb)\, and the associated 
distributions P (A t ) and P(A m ). We monitored 
these functions after deducting a trivial noise |llj . 
namely, small thermal events of typical size |4()| : 



A = 



ln[l+p+(l-p)e- 2 \ h \/ T ]' 



(2) 



These events correspond to the short-length fluctu- 
ations that are evident by direct inspection in the 
inset of figure 1 . 

The distribution P (A m ) that results after de- 
ducting small events is illustrated in figure 2. This 
nicely fits 

P(A ro )~A--W (3) 
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Figure 2: Log- log plot of the size distribution 
P (A m ) of large avalanches for an ensemble of in- 
dependent particles of radius (from bottom to top) 
R = 30, 42, 60, 84 and 120, respectively. Plots of 
the duration distribution P(A t ) vs. cA] for each 
R are also shown (x), with c w 0.5 and 7 « 1.52 
(see text). For visual convenience, the curves are 
shifted vertically by 4™ with n = to 4 from bot- 
tom to top. Running averages have been performed 
for clarity purposes. 



with 



r(R) 



a\R 



(4) 



where Too = 1.71 (4) . Figure 3 depicts the corre- 
sponding duration distributions. They follow 



-«CH) 



with 



a (R) = a 

.DC + a 2 R 



(5) 



(6) 



where a OT = 2.25 (3) . In both cases, size and du- 
ration, the apparent power law ends with an expo- 
nential tail, 



P(A) ~ exp(-A/A* 



(7) 



The cutoffs that we observe follow A* ~ R@ with 
P m ~ 2.32(6) and (3 t ~ 1.53(3), respectively (see 
inset to Fig. 3). 

We also determined that observing a power law 
requires both free borders and the noncquilibrium 
condition. That is, the distributions P (A) look ap- 
proximately exponential if the system has periodic 
borders and/or one sets p — in eq. QJ. This is 



Figure 3: Log- log plot of the duration distribution 
P (Aj) for the same ensembles of particles as in fig- 
ure 2. For visual convenience, the curves are shifted 
vertically by 2™ with n — to 4 from bottom to top. 
Running averages have been performed for clarity 
purposes. Inset: log-log plot of the size (top) and 
duration (bottom) cutoffs A* vs. R. Lines are 
power-law fits. 



discussed below. Another main result is that the 
observed apparent power laws are here a sum of 
exponential contributions. 

To prove the latter result, we followed the de- 
magnetization process in a large circular particle. 
The main interest was in the interface between the 
rich and poor spin-up regions at low temperature. 
One observes curved interfaces due to the faster 
growth of the domain near the concave open bor- 
ders. In fact, the critical droplet always sprouts at 
the free border [41| . Then, given that curvature 
costs energy, the large avalanches tend to occur at 
the curved walls, which then transform into rather 
flat interfaces. We confirmed this by estimating the 
mean avalanche size (A m ) and interface curvature 
(C) as a function of magnetization m. The curva- 
ture C is defined here as the number of rising steps 
at the stable-metastable interface We plot in 
figure 4 our results for these observables. After av- 
eraging over many runs, definite correlations show 
up. That is, as one could perhaps have imagined, 
the event size is determined by the interface curva- 
ture just before the avalanche occurs. 

This is confirmed by monitoring P (A m | C) , 
the conditional probability that an avalanche of 
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Figure 4: Semilogarithmic plot of (A TO ) (solid line) 
and (C) (dotted line) as a function of magnetiza- 
tion, after averaging over 3500 independent runs. 
Notice the non-trivial structure uncovering a high 
degree of correlation between the mean size of 
avalanches, (A m ), and the average curvature (C) 
of the interface at which the avalanche originates. 



Figure 5: Semilogarithmic plot of P (A m | C) , the 
size distribution for avalanches developing at a wall 
of constant curvature, C; C increases from bottom 
to top. Here, N m = ^NA m . (For visual conve- 
nience, the curves are shifted vertically by 10™ with 
n = to 4 from bottom to top.) Running averages 
have been performed for clarity purposes. 



size A m develops at an interface region of curva- 
ture C. We studied this in great detail by simu- 
lating an interface of constant curvature evolving 
by Q ■ Figure 5 shows that P ( A m | C) has two 
regimes for given C. The first one corresponds to 
the small thermal events mentioned above, namely, 
those of typical size given by (0). The second 
regime exhibits, contrary to the situation in fig- 
ure 2, (stretched ) exponential behavior, namely 
P(A m | C) ~ exp [- (A m /A m ) n ] with 77 a 0.89. 
That is, a wall of curvature C induces avalanches 
of typical size A m (C) . 

This fact turns out most relevant because, due 
to competition between the randomness induced by 
free borders and the one induced by p in (JJJ, the 
interface tends to exhibit a broad range of curva- 
tures with time. More specifically, relaxation pro- 
ceeds via a series of different configurations, each 
characterized by a typical curvature of the interface 
and by the consequent typical form of the critical 
droplet inducing the avalanche. Therefore, what 
one really observes when averaging over time is a 
random combination of many different avalanches, 
each with its typical well-defined (gap-separated) 
size and duration, which results in an effective dis- 
tribution. The fact that this combination depicts 
several decades (more the larger the system is) of 



power-law behavior can be understood on simple 
grounds. 

Let Q (A) the probability of A, and P (x \ A) = 
Aexp(— Ax) the probability of an event of size x 
given A. Assume that A can take a finite number 
of equally spaced values A^, k — 0, 1, 2, . . . , n, in 
the interval [A m { n , ^4 max ] , so that A^ = A mnl + kd 
with 5 = (j4 ma x — ^4min) /n (alternatively, one may 
assume randomly distributed A^s), and that all of 
them have the same probability, Q (A) = const. 
One obtains that 



P(x) 



Se -xA min f 



1 - e 



-x5 



A T 



Ar, 



-(n+l)xS 



-nxd ' 



1 



(8) 



The fact that even such a simple, uncorrelated 
ansatz describes qualitatively the data is illustrated 
in figure 6. That is, the superposition of a large 
but finite number of exponential distributions, each 
with a typical scale, yields an effective global distri- 
bution which is consistent with apparent scale in- 
variance. This distribution extends in practice up 
to a cutoff, which is also observed in experiments. 
This cutoff, which corresponds to the slowest expo- 
nential relaxation, is given by exp (— xA m { n ). There 
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is no evidence that a more involved computation 
would modify this qualitative conclusions. How- 
ever, taking into account dynamic correlations as 
revealed by figure 0] is certainly needed in order 
to improve quantitative predictions. In particu- 
lar, eq. (jHJ predicts a size-independent exponent 
t(R) = Too = 2, somewhat different from the ob- 
served asymptotic t m = 1.71(4). 

Consider now P(A t | A m ), i.e., the probabil- 
ity that the avalanche of size A m lasts a time A t . 
We confirmed that this exhibits well-defined peaks 
corresponding to large correlations, i.e., avalanches 
of a given size have a preferred duration and vice 
versa. Assuming A m ~ A], we obtain 7 
I,,: . >■ — 1.52(5). Using this relation, one may 
obtain the duration distribution by combining I© 
with P(A m )dA m = P(A t )dA t . A comparison of 
the resulting curve with data in Fig. H3 leads to 
7 ~ 1.52, in agreement with the value obtained 
from the cutoff exponents (3. More generally, a 
scaling plot of P(A t ) vs. cA] , with c some pro- 
portionality constant, must collapse onto the cor- 
responding curve P(A m ) for each R. This is con- 
firmed in Fig. 2 for 7 ~ 1.52, further supporting 
the scale-superposition scenario. On the hypothesis 
that both P(A m ) and P(A t ) were true power-law 
distributions, one would obtain the scaling relation 
(a — 1) = j(t — 1). However, our values above for a 
and t would imply here that 7 ~ 1.76. This misfit 
is a consequence of the fact that, according to our 
point in this paper, none of the distributions P(A) 
exhibits true scaling behavior. 

4 Discussion 

The model studied in this paper has been demon- 
strated to relax only via well-defined near- 
exponential events. Each has its own scale, but 
many of them randomly combine into a distribu- 
tion that exhibits a power-law portion. This occurs 
in the model because two of its features, namely, 
p =/= and free borders. Otherwise, i.e., p = 
and/or periodic boundary conditions, the appar- 
ent power law does not emerge, but one observes a 
well-defined mean. Therefore, the cause for a rela- 
tively broad range of possible different scales is the 
underlying nonequilibrium randomness that char- 
acterizes the model. The question is whether this 
picture applies also to reported scale-free fluctua- 




Figure 6: Solid lines are predictions from equation 
© for n = 200, A min = 0.007 and A max = 1. The 
symbols stand for the avalanches duration (lower 
curve) and size (upper curve) when R = 60, i.e., 
two of the data sets in figures 2 and 3. In this 
particular case, the finite-size exponent is t(R = 
60) = 2.06(2), allowing direct comparison with eq. 

©■ 



tions in many natural phenomena, where it is diffi- 
cult to investigate separate elementary events. We 
show below that there are some indications — but 
not a proof — that this may be the case. We also 
remark in this section some important features of 
our picture that one should look for in experiments. 

We first remark that analyzing the fluctuations 
of interest requires a previous scrutiny of data in 
the model separating small from larger events. The 
former are random events of well-defined size ac- 
cording to eq. @ ■ The latter are more structured, 
correlated events or "avalanches" . This separation 
is theoretically motivated > an d it is also sup- 
ported by experiments |28l I43j . The avalanches are 
then accurately described by an apparent power 
law. This involves an exponent and a cutoff both 
depending on the size of the system. Interesting 
enough, there is a well-defined limit for a macro- 
scopic particle. That is, even though free borders 
(a surface effect) are essential to the phenomena, 
increasing the particle to macroscopic will not pre- 
vent one from observing all scales, including very 
large, say macroscopic avalanches. 

Failing in separating accurately small from large 
events will, in general, result in a non-significative 
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distribution. On the other hand, though the de- 
scribed behavior is robust within a broad region 
of parameter space, unambiguously observing the 
relevant phenomena requires some care in order 
to have compact enough clusters and well-defined 
interfaces. As indicated in section ^ prevalence 
of one ingredient over the others may importantly 
hamper statistics which may obscure the situation. 
These two facts — together with other agents EU 
(see below) — seem to produce both power-law and 
exponential distributions in closely related situa- 
tions as described, for instance, in [2]. 

There is no indication that the apparent scale 
invariance that we observe is to be associated with 
chaos, e.g., sensitivity to the initial condition. Our 
picture is neither consistent with a critical con- 
dition. Criticality implies a diverging correlation 
length which is not detected in our simulations. 
Instead, the many scales in the model are sim- 
ply due to the underlying nonequilibrium random- 
ness mentioned above. A balance between ran- 
dom and nonequilibrium features has already been 
claimed to be essential for SOC in different set- 
tings CHI ES EH EH EH Our study suggests 
this is the essential physics in a family of situa- 
tions. That is, even though details may vary, e.g., 
from random interface rearrangements in magnets 
to slip complexity in earthquakes, there is always 
some microscopic randomness which induces mul- 
tiple short-lived situations. This constantly halts 
the decay and, ultimately, leads to apparent scale 
invariance. 

We find it remarkable — even though it does 
not prove anything — that the statistical proper- 
ties of the avalanches in our model are indistin- 
guishable in practice from what has been reported 
for some detailed laboratory experiments. For in- 
stance, size corrections similar to the ones in (J2J 
and © for r and a, respectively, have been re- 
ported in avalanche experiments on rice piles |44j . 
Moreover, our values for the infinite case are strik- 
ingly close to the ones reported in magnetic ex- 
periments, e.g., ^ — 1-77(9), = 2.22(8) and 
7 = 1.51(1) in Ref. |11| for quasi-two dimensional 
systems. (See also Ref. E2|.) O n the other hand, 
our cutoff values in @ follow the precise trend ob- 
served in magnetic materials 0B1EH]- More quali- 
tatively, one may argue that the available literature 
on the Barkhausen noise, which is rather proclive to 
a hypothetical critical point, provides meaningful 



indications of consistency with our non-critical be- 
havior, as reported elsewhere EH- If this is granted, 
the many transitions through halt states in our 
picture would correspond to topological rearrange- 
ments of domain walls in the Barkhausen case. In 
fact, there is also in this case a domain finite size 
and a nonequilibrium perturbing drive (the varying 
field) , that would induce a condition similar to our 
"nonequilibrium randomness" constantly changing 
the scale. 

Avalanches in our model do not continue in- 
definitely in time but disappear when the system 
reaches the stable state. This contrasts in principle 
with the stationary character of some experimental 
signals. However, it does not prevent our compar- 
ison. To illustrate this, imagine that we flip the 
magnetic field sign in our model every time the sys- 
tem reaches the final stable state. This would give 
rise to a cyclic steady process, much in the spirit 
of hysteresis experiments in magnets. Avalanches 
observed in this cyclic state have the same proper- 
ties that the ones reported here. In particular, our 
conclusions on the origin of the apparent scaling of 
distributions remain unchanged, while the process 
is now cyclically stationary. 

No doubt it would be interesting to study the 
possible occurrence of "short-lived halt states" in 
nature. These are associated in the model with flat 
interfaces. That is, once the initial metastability 
breaks down, the particle becomes inhomogeneous, 
and flat interfaces have a significant probability to 
form after each avalanche (which aims at minimiz- 
ing interfacial energy). As this is the most stable 
configuration against small perturbations, the sys- 
tem remains some time with constant magnetiza- 
tion m (tb) ■ This may be described as an entropic 
metastability. There is no real energy barrier but 
an unstable situation such that a given microscopic 
random event suffices to initiate the next avalanche. 

Interesting enough, this picture gives more hope 
to the goal of predicting large events. That is, 
we claim that catastrophes are not a rare random 
emergence in an strongly correlated bulk which, 
consequently, have the same cause as the small 
events. |2 El E] Instead, events are characterized 
by their size, and each size follows from some spe- 
cific microscopic configuration. The configurations 
that, under appropriate conditions, may originate 
large events qualitatively differ from the ones cor- 
responding to smaller events. In summary, there is 
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some specific cause for each event which depends 
on its size. Consider for instance a stock market; 
its evolution is also characterized by discontinu- 
ous, sudden jumps between different locally-stable 
states. In order to predict a crash in this system, as- 
suming that our picture applies, one should look for 
the simultaneous occurrence of a large flat interface 
(predisposition of the players) and macroscopic free 
borders (some large external perturbation) . Study- 
ing the statistical properties of the many events will 
then only inform on the relative probability of each 
microscopic relevant configuration. 

A detailed description of the scape mechanism 
from entropic metastability could relate our picture 
to other approaches. We mention in this respect 
that, after averaging over many independent par- 
ticles, the lifetimes of the halt states during the 
relaxation of our model depict an exponential dis- 
tribution. Therefore, they show a typical scale. 
This scale turns out to be much shorter than the 
time scale for the system relaxation, as reported to 
characterize the supposed nonequilibrium critical- 
ity which is assumed to under ly many 1// noises; 
see, for instance, Ref.0- An even more detailed 
look, which requires averaging over time intervals, 
reveals in our model that this scale definitely de- 
creases with t from, say, macroscopic (> 10 5 MCS) 
to microscopic (~ 10 MCS). This feature, which is 
already evidenced by (direct inspection of) figure 
1, is one that could perhaps be easily detected in 
experiments. 

Finally, we remark that there are other possi- 
ble explanations for 1// noise based on non-critical 
mechanisms; see, for example, Refs.|21 ) -|26 | . These 
are less general than the mechanism proposed here, 
and often restricted to some very specific situation. 
Furthermore, some of these descriptions may be in- 
terpreted at the light of a superposition of many dif- 
ferent typical scales, as in our mechanism. A simi- 
lar origin for electronic 1// noises was suggested in 
the past (see, for instance, R.efs. [T7l 115111911201 l2"T] V 
though the present paper is, to our knowledge, 
the first one in establishing an explicit relation be- 
tween elementary events (avalanches) and micro- 
scopic physical processes. 
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